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Abstract 

The quasi - relativistic harmonic oscillator bound - states con- 
structed by Znojil ( 1996 J. Phys. A29 2905) are investigated via a 
new methodical proposal. Compared to those obtained by an anony- 
mous referee ( from a direct numerical integration method) of ZnojiPs 
paper [3], our results appeared to be more favorable than those ob- 
tained by Znojil via quasi - perturbative, variational, Hill - determi- 
nant and Riccati - Pade methods. Bound - states with larger angular 
momenta I are also constructed. 



Among exactly soluble Hamiltonians exists the harmonic oscillator (HO) 
Hamiltonian 



H (HO) = P_ + l _ mu 2 Y 2_ (1) 

2m 2 K J 



The equidistant form of its spectrum have attracted attention in quantum 
control theory [1] and represents, in view of [2], a ' long sought - after dream' 
of ' steering wavepackets into desired states '. Moreover, it fits, via exper- 
imental observations, vibrational excitations of molecules and some low - 
lying energy levels in atomic nuclei. 

On the other hand, a free spin - field fulfills Klein - Gordon equation 



(E 2 - H 2 )^ = ; H = \/c 2 p 2 + m 2 c 4 , (2) 

or the Schrodinger formulation of it (E — H)^/ = 0. Of course there exist two 
admissible solutions for a given momentum p, i.e. positive and/or negative 
energy solutions. However, we devote our work to the positive solution in 
the Schrodinger - formulated Klein - Gordon equation. With the minimal 
coupling, a Lorentz 4 - vector harmonic oscillator potential is coupled as 
the - component of the four - vector potential, i.e. E — >■ E — eA with 
eA = mcu 2 r 2 /2. Hence, a fully relativistic description of H^ H °^ leads to the 
emergence of a quasi - relativistic harmonic oscillator (QHO) Hamiltonian 

H (QHO) = ^J m 2 c 4 + p 2 c 2 + I mCc ,2 r 2. (3) 

Which, in momentum representation, implies the one - dimensional Schrodinger 
equation 



—mh 2 uj 2 
2 



dp 2 



1(1 + 1) 
p 2 



^(p) + [Jm 2 c 4 + p 2 c 2 -mc 2 }^(p) =EV(p), (4) 



where r 2 = — h 2 A p , p G (0, oo), E is the binding energy and / = 0, 1, • • • 
denotes the angular momentum quantum number. Of course, on the asymp- 
totically physical grounds, the wavefunction ^(p) satisfies the boundary con- 
ditions [3] 



*(p) 



pi+i j or |pj << i 
exp(—vp 3 / 2 ) , v - 



:(hw) \ for \p\ » 1. 



A rescal of the variable p through p = ymKuq would, in turn, lead to a 
transparent form of Schrodinger equation [3] . Strictly, 



1 d 2 1(1 + 1) 



+ 



2 dq 2 2q 



V(q) 



^n r ,l(0) = £ n r ,l^n r ,l(<l), 



(5) 



where, e nr>l = 2E nrtl /(hu), 



V(q) 



Q' 



:(V S RAo(q) ~ 1), 



(6) 



VsRAo(q) = Jl + a 2 q 2 , 



(7) 



a 2 = huj/(mc 2 ), and the square - root anharmonic oscillator potential (7) 
simulates a quasi - relativistic squeezing of the harmonic oscillator spectrum. 
The transition from the parabolic HO well, in standard coordinate repre- 
sentation, to the hyperbolic shape Vsrao(q)i m momentum representation, 



is claimed to prove phenomenologically useful and methodically challeng- 
ing. Znojil [3] has, therefore, invested several eligible ( namely, perturbative, 
variational, Hill - determinant, and Riccati - Pade) methods to construct its 
bound - states. With the permission of an anonymous referee of his paper, 
Znojil has also reported ( in table 2 of [3]) the referee's results from direct 
numerical integrations. 

Whilst using a quasi - perturbation prescription ( Eq.(ll)in [3]), a loss 
of precision occurred at "large" a — 1/2 ( table 1(b) in [3]). Upon which 
the anonymous referee remarked that it may also cause a loss of the upper - 
bound character of the quasi - perturbation prescription. In accordance with 
a second referee's remark, being curable by a Pade - type resummation, the 
loss of boundedness phenomenon may emerge at any a. 

To the best of our knowledge, the paper of Znojil [3] is the only available 
one in the literature, for the quasi - relativistic harmonic oscillator, and 
merits therefore further considerations. 

In this paper we formulate a new method to solve the Fourier transformed 
Schrodinger equation (4), with V(q) represented by (6). Our method consists 
of using l/l as an expansion parameter, where I — I — /3, I is a quantum 
number, and (3 is a suitable shift introduced, mainly, to avoid the trivial 
case I = 0. The spiritual soundness of the textbook perturbation theory is 
therefore engaged. Hence the method should be called pseudoperturbative 
shifted - / expansion technique (PSLET). 

With the noninteger ( irrational) orbital angular momentum I, equation 
(5) reads 

{Id 2 ~ 1 

1 2 df + V{ - q) j * nr ' 1 ^ = £n r^n„M, (8) 



v-fa) = P±M±g±W±i) + P , /( „. (9) 



Where Q is a constant that scales the potential V(q) at large - / limit and is 
set, for any specific choice of / and n r , equal to I 2 at the end of the calculations 
[4-9]. And, j3 is to be determined in the sequel. 

Our systematic procedure begins with shifting the origin of the coordi- 
nate through 

x = l 1 /\q-q )/q , (10) 

where q is currently an arbitrary point to perform Taylor expansions about, 
with its particular value to be determined. Expansions about this point yield 



2_ = f*(_ 1) n(fi+ll x nj-n,2 i 



o-z^ ■ > •• ■ (H) 



= g r^m\ (^r r _„ /2 _ (12) 



V(x 



It should be mentioned here that the scaled coordinate, equation (10), has no 
effect on the energy eigenvalues, which are coordinate - independent. It just 
facilitates the calculations of both the energy eigenvalues and eigenf unctions. 
It is also convenient to expand e nr j as 

oo 

e nr ,i = E ^r n - (13) 

n=-2 



Equation (8) thus becomes 



\& + T*«M 



q 2 



(14) 



with 



I 



V(x(q)) = q 2 J 



1 + V(q ) 



2g 2 Q 



+ l 1 ' 2 



-x + 



Q 



+ 



AJ2 



-x z + 



+ (2/3+i)E(-i; 



+ ?o 2 E 



n=3 



V (Qo)Q qX 
2Q 

' (n+22 

, 2g D 2 ^V dQo J n\Q 



(n + 1) 



-x 



ni—n/2 



71=1 

cH/(g )\ (g x) n 



J-(n-2)/2 



+ w+ i)g(_i)»(!L^i a; »r-(»+2)/2 + (^+l) 



(15) 



n=0 



where the prime of V(g ) denotes derivative with respect to q . Equation 
(14) is exactly of the type of Schrodinger equation for one - dimensional an- 
harmonic oscillator 



W + ^ V + A » + PW 



X nr (x) = \n r X nr (x), 



(16) 



where P(x) is a perturbation - like term and A is a constant. A simple 
comparison between Eqs.(14), (15) and (16) implies 



An = l 



2 Q 



, 2/3+1 /?(/?+!) 
+ 2 + 2^ 



(17) 



2 Q 



+ 



2/3 + 1 , L_ 

-^— + (n r + -)n 



+ 



W±il + A <o, 

2 r 



+ e ASr 1 ^-", 



n=2 



(18) 



and 



a - o 2 V f (n) H n+1 ) 



n=-2 



Equations (18) and (19) yield 



(19) 



-(-2) 

"n r ,l 



i + n<? ) 



2g Q 2 g 



(20) 



-(-1) = ±_ 

"n r ,l 2 



-^— + (n r + -)fi 



(21) 



jo) = J_ 

n r ,l 2 



M±l) +A (o) 



(22) 



© = *$?/<£ ; n > l. 



(23) 



Here q is chosen to minimize e n /, i. e. 



—^- = and ":■' > 0, 

dq Q dq 2 



(24) 



which in turn gives, with / = ^JQ 



l-P = ^q 3 V'(q ). 



(25) 



Consequently, the second term in Eq.(15) vanishes and the first term adds a 
constant to the energy eigenvalues. 

The next leading correction to the energy series, le~ l , consists of a con- 
stant term and the exact eigenvalues of the unperturbed harmonic oscilla- 
tor potential Q 2 x 2 /2. The shifting parameter f3 is determined by choosing 
le~ l — 0. Hence 



P 



- + (n r + -)Q 



(26) 



where 



n 



3 + 



QoV"(q ) 



\ V'(q ) 



(27) 



Then equation (15) reduces to 



jV(x(q)) = q 2 l 



1 V(qo) 
2q 2 o Q 



( n )(^M~ n / 2 



n=0 



(28) 



where 



V l°)( x ) = ±tfa* + W±±, 



(29) 



V W(x) = -(2/3 + l)x - 2x 3 + q ° V (go) x 3 , 

QQ 



(30) 



and for n > 2 



v {n) (x) = (- 



4 )"(2/3 + 1)^±^U" + (-l)" ^ +1) (n - l)a>" 2 ) 



+ 



1 J 2 Q(n + 2)! dg^° 



-.n+2 



X 



n+2 



(31) 



Equation (14) thus becomes 



1 J2 oo 

£ ux n=0 

1/ W+l) 
/ I 2 



*«_.i(a:) 



« + e Air 1 *-* 



n=2 



^n r ,/(a:) 



(32) 



When setting the nodeless, n r = 0, wave functions as 



Woi(x(q)) = exp(U ,i(x)), 



(33) 



equation (32) is readily transformed into the following Riccati equation: 



-1 oo 

- -\U'\x) + U\x)U\x)} + Y, v {n \x)l- n/2 

^ n=0 



1 f 13(13 + 1) 



I 



+ ^o 



(0) 



£ Xr L) r n - (34) 



n=2 



Hereinafter, we shall use U(x) instead of Uqj(x) for simplicity, and the prime 
of U(x) denotes derivative with respect to x. It is evident that this equation 
admits solution of the form 



10 



u'( x ) = E u {n \x)r n / 2 + E G {n \x)r^ 



l)/2 



n=0 



n=0 



(35) 



where 



n+l 



C/W( x ) = ^ £>. 



.2m-l 



m,n<K i Uo,n U, 



m=0 



n+l 



G W W = E Cm,r 



x 



2m 



m=0 



Substituting equations (35) - (37) into equation (34) implies 



(36) 



(37) 



1 oo 

-E 

^ n=0 



£y(n) j-n/2 _j_ gr(n) j-(n+l)/2 

I oo oo 
_ V" V" \jj(n)jj(p)J-(n+p)/2 _|_ Qf(n)g f (p)^-(n+p+2)/2 + 2 [/(n) Qf(p)p(n+p+l)/2 

(36 



n=0p=0 

/ I 2 



+ e ^" n/2 = t ( ^r^ + A o 0) + E Ar i} r n , 



n=0 



n=2 



where primes of C/^(x) and G^^rr) denote derivatives with respect to x. 
Equating the coefficients of the same powers of I and x, respectively, ( of 
course the other way around would work equally well) one obtains 



_Itf(o)'_Ij7(o)i7(o) +v (o) =0 
2 2 



(39) 



U^\x) = D lt0 ; D lfi = -Q, 



(40) 
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and integration over dx yields 

U i0 \x) = -Six. (41) 

Similarly, 

_ I[[/(i)' + G (o)'] _ 17(0)17(1) _ £7(0)^(0) + v (i) = 0> (42) 

C/ (1) (x) = 0, (43) 

G(°)(a;) = Co,o + C 1 , x 2 , (44) 



^ - -£. (^) 



C 0) o = ^(C7 li o + 2/3 + 1), (46) 



— J§^ <-' 



--[c/( 2 )' + g«'] --J2 c/Hc/( 2 -) - Ig(°)g(°) 

2 2 n=0 2 
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- J2 U^G^ + v^ = /3(/3+1) + \$\ (48) 

n=0 ^ 



U^ 2 \x) = D h2 x + D 2 ^x\ (49) 



G (1) (a;) = 0, (50) 



1,<3 



2 



,0 



^2,2 = ^(^-^2) (51) 



£>i,2 = ^(^2,2 + C 0) oCi )0 -|(2/?+ 1)), (52) 



5 g p 6 cfV( go ) 
2 + 24Q dqt 



B. = ^^~^, (53) 



Ai 0) = -2(^ + Co 2 o)- (54) 

• • • and so on. Thus, one can calculate the energy eigenvalue and the eigen- 
functions from the knowledge of C m ^ n and D mn in a hierarchical manner. 
Nevertheless, the procedure just described is suitable for systematic calcula- 
tions using software packages (such as MATHEMATICA, MAPLE, or RE- 
DUCE) to determine the energy eigenvalue and eigenfunction corrections up 
to any order of the pseudoperturbation series. 

13 



It should be mentioned that the energy series, Eq.(13), could appear 
convergent, divergent, or asymptotic. However, one can still calculate the 
eigenenergies to a very good accuracy by forming the sophisticated Pade ap- 
proximants to the energy series [11]. The energy series, Eq.(13), is calculated 

u p to 4? A 8 b y 



i2„(-2) 



-(0) 



(8)/78 



ev = l'etf } +eW + ~- + e%/r + 0(l/n 



(55) 



and with the P$(l/l) and P N (1/0 Pade approximants it becomes 



72(-2) . jjN/ 



eoAN,N]=l'e)TS ) +mi/l). 



(56) 



and 



eoAN,N+l]=Pei ) - l 2) +P^ 1 (l/l). 



(57) 



Our strategy and prescription are therefore clear. 

Let us now consider the Fourier transformed equation, Eq.(4), with the 
rescaled variable p = y/mhuq, represented by equations (5)-(7). The substi- 
tution of Eq.(6) in (26), for n r = 0, implies 



(3 = --(i + n) ; n 



4 + 3a 2 q* 
\TT^qJ 



(58) 



Eq.(25) thus reads 



\ 



1 + a 2 q, 



2 nZ 



(59) 
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Equation (59) is explicit in q and evidently a closed form solution for q is 
hard to find, though almost impossible. However, numerical solutions are 
feasible. Once q is determined the coefficients C myn and D myTl are obtained 
in a sequential manner. Consequently, the eigenvalues, Eq.(55), and eigen- 
functions, Eqs.(35)-(37), are calculated in the same batch for each value of 
a, and /. 

In order to make remediable analysis of our results we have calculated the 
first ten terms of the energy series. The effect of each term has been taken 
into account. We have also computed the Pade approximants Eoj[N,M] for 
N = 2,3,4 and M = 2,3,4,5. Therefore, the stability of the energy series 
and that of the sequence of Pade approximants are in point. 

Table 1 shows PSLET results for the ground - state energies £o,o> covering 
a wide range of the anharmonicity a, along with the exact ones from direct 
numerical integration method ( DNI), carried out by the anonymous referee 
of paper [3]. To avoid exhaustive numbers of tables we do not list Znojil's 
results. However, we will just refer to them. A comparison between PSLET 
and DNI results implies excellent agreements. The nice trend of stability of 
the energy series (55) (i.e., a signal of nice course of convergence.) is well 
pronounced. The effect of the higher - order corrections on the first few terms 
of the energy series bears this out. 

In contrast with Znojil's results ( table 1(b) in [3]) for "large" a=l/3 and 
a=l/2, via a quasi - perturbative prescription ( Eq.(ll) in [3]), there is no 
indication that our series will blow up at higher - orders and ( since, mainly) 
our expansion parameter (1//) is less than one for all values of a reported in 
the text. Of course there is always the contribution of the (K+ 1) - term, but 
so far our prescription is performing so good. Whilst Znojil's prescription 
marks nice stability for small a, sever oscillations of his series occur at low 
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- order, especially for "large" a— 1/2, causing in effect a breakdown in the 
boundedness character of his prescription. Although this phenomenon is 
curable by resummation tools like the sophisticated Pade approximants, as 
suggested by the second referee of his paper [3], however, this, in our opinion, 
shall not dramatically cure the loss of precision in Znojil's results ( table 1(b) 
in [3]), to be documented in the sequel. 

Switching to alternative methods for independent checks of his numerical 
predictions, Znojil used the Hill determinant and Riccati - Pade methods. In 
the light of his experience in the Hill determinant, an onset of convergence is 
clearly manifested ( table 3 in [3]), but larger dimensions and/or improved 
elementary convergence factor would be necessary to reach the domain of 
more satisfactory numerical precision. Upon his experience, moreover, in a 
slightly more complicated ( compared to the Hill determinant) Riccati - Pade 
method ( RPM), the 11- dimensional Toplitz determinants have offered very 
satisfactory precision ( table 4 in [3]). However, a typical bizarre character- 
istic of the RPM is well documented [3,11]. Namely, it leads to a number 
of clustered solutions, for a given value of the coupling a, resulting from the 
existence of several eligible physical roots of the Hankel [11] or Toplitz [3] de- 
terminants. Yet the ambiguity of these roots increases with the dimensional 
growth of the determinants. Although clustering is a good signal of being 
close to a physical root, a decision on which of these roots is the best has to 
be made. So far, to the best of our knowledge, a general way of establishing 
this property has not been found. 

The effect of the angular momentum quantum number / on the stability, 
hence on convergence, and precision is reported in table 2 for a=l/2. Con- 
fidently, one concludes that better convergence and more precise numerical 
results are obtained as / grows up. Similar effect should be expected from 
the nodal quantum number n r ; I and n r have almost identical effects on our 
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pseudoperturbative expansion parameter I. 

The stability of the sequence of Pade approximants ( table 3) is fascinat- 
ing. Although there is no signal that our series will blow up, table 1 marks 
this fact, the effect of Pade approximants on precision is limited. For a fixed 
a, say 1/2, more precision is obtained via Pade approximants as I increases 
( table 4). Adhered to the conventional practice of perturbative calculations 
( i.e. only a few terms of a " most useful" perturbation series reveal the 
important features of the solution before a state of exhaustion is reached.), 
we list PSLET results ( table 5) from the first 6 terms of our energy series 
with £o,o[3, 3] Pade approximant. Compared to those from DNI, our results 
are readily satisfactory. 

To summarize, we have used a new methodical proposal to investigate 
the bound - states of the quasi - relativistic harmonic oscillator. Using the 
perturbation expansion parameter 1/7, we have demonstrated that our ap- 
parently artificial perturbation recipe PSLET is convincingly powerful and 
methodically practical. 

Perhaps it should be noted that for each entry in tables 1-5 one can 
construct the wavefunction from the knowledge of C myn and D myTl . However, 
such a study lies beyond the scope of our methodical proposal. 

In addition to Znojil's interpretation of the QHO Hamiltonian (3) that it 
leads, in effect, to the formally correct relativistic Dirac equation, we have 
shown that it could as well represent a Klein - Gordon particle in a parabolic 
well. Precisely, the 4 - vector potential eA = —Ze 2 (3 — r 2 /a)/2a, or in 
short eA = A + Br 2 , represents an improved approximation for a realistic 
pionic atom, hence the Hamiltonian in Eq.(3) addresses the Klein - Gordon 
Hamiltonian for the potential of a homogeneously charged sphere [13]. 

The applicability of our recipe extends beyond the present quasi - rel- 
ativistic harmonic oscillator model. Some applications are in order. The 
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eigenstates of a hydrogenic impurity in a spherical quantum dot (QD) [14]. 
Quasi - two - dimensional QD helium [15]. Two - electrons QD in a magnetic 
field [16]. Excitons in harmonic QD [17]. Hydrogenic impurity or heavy 
excitons in arbitrary magnetic field [18,19], • • • etc. 
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Table 1: Ground - state energies , where K represents the first K - terms 
of Eq.(54) and DNI from direct numerical integration [3]. 



K 


a = 1/100 


a = 1/20 


a - 1/10 


a = 1/5 


1 


2.999906 257850 8 


2.997661 14471 


2.990702 721 


2.963706 98 


2 


2.999906 259959 4 


2.997662 45159 


2.990723 092 


2.964001 14 


3 


2.999906 259959 1 


2.997662 44641 


2.990722 775 


2.963984 17 


4 


2.999906 259959 1 


2.997662 44644 


2.990722 783 


2.963985 58 


5 


2.999906 259959 1 


2.997662 44644 


2.990722 782 


2.963985 42 


6 


2.999906 259959 1 


2.997662 44644 


2.990722 782 


2.963985 45 


7 


2.999906 259959 1 


2.997662 44631 


2.990722 775 


2.963985 06 


8 


2.999906 259959 1 


2.997662 44635 


2.990722 777 


2.963985 17 


9 


2.999906 259959 1 


2.997662 44635 


2.990722 777 


2.963985 17 


10 


2.999906 259959 1 


2.997662 44634 


2.990722 776 


2.963985 14 


DNI 





2.997662 44644 


2.990722 782 


2.963985 44 




a- 1/4 


a = 1/3 


a = 1/2 


a = 2 


1 


2.944289 62 


2.904543 


2.80482 


1.9189 


2 


2.944955 82 


2.906342 


2.81083 


1.9331 


3 


2.944899 04 


2.906100 


2.80951 


1.9334 


4 


2.944905 92 


2.906145 


2.80987 


1.9323 


5 


2.944904 81 


2.906135 


2.80977 


1.9319 


6 


2.944905 03 


2.906138 


2.80979 


1.9322 


7 


2.944903 79 


2.906132 


2.80978 


1.9328 


8 


2.944904 11 


2.906134 


2.80977 


1.9329 


9 


2.944904 11 


2.906134 


2.80977 


1.9329 


10 


2.944904 03 


2.906133 


2.80975 


1.9315 


DNI 


2.944904 99 


2.906136 


2.809786 


1.932334 
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Table 2: The effect of the angular momentum quantum number I on con- 
vergence and precision for a = 1/2. 



K 


/ = 1 


1 = 5 


1 = 10 


1 = 20 


1 


4.569573 


10.977406 1 


17.963301 722 


29.948955 412735 


2 


4.575998 


10.981520 1 


17.965505 044 


29.949880 189612 


3 


4.575031 


10.981360 3 


17.965487 822 


29.949881 669238 


4 


4.575173 


10.981349 8 


17.965484 484 


29.949881 347421 


5 


4.575167 


10.981352 9 


17.965484 566 


29.949881 340115 


6 


4.575155 


10.981352 8 


17.965484 596 


29.949881 340512 


7 


4.575164 


10.981352 7 


17.965484 595 


29.949881 340559 


8 


4.575161 


10.981352 7 


17.965484 595 


29.949881 340559 


9 


4.575161 


10.981352 7 


17.965484 595 


29.949881 340559 


10 


4.575162 


10.981352 7 


17.965484 595 


29.949881 340559 



Table 3: The effect of Pade approximants on convergence and precision. 



£o,o [N,M] 


a = 1/100 


a= 1/10 


a = 1/3 


a = 2 


£o,o[2,2] 


2.999906 259959 1 


2.990722 78 


2.90614 


1.9319 


£o,o [2, 3] 


2.999906 259959 1 


2.990722 78 


2.90614 


1.9323 


£o,o[3, 3] 


2.999906 259959 1 


2.990722 78 


2.90614 


1.9334 


£o,o[3,4] 


2.999906 259959 1 


2.990722 78 


2.90614 


1.9325 


£o,o[4,4] 


2.999906 259959 1 


2.990722 78 


2.90613 


1.9326 


£o,o[4,5] 


2.999906 259959 1 


2.990722 78 


2.90613 


1.9326 
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Table 4: The effect of / and Pade approximants on convergence and precision 

for a = 1/2. 



£o,i[N,M] 


/ = 


1 = 1 


1 = 3 


eo,i[2,2] 


2.809788 


4.575157 


7.893680 


eo,i[2,3] 


2.809826 


4.575161 


7.893681 


£ cy[3, 3] 


2.809783 


4.575161 


7.893681 


£o,j[3,4] 


2.809767 


4.575161 


7.893681 


eo,i[4,4] 


2.809790 


4.575161 


7.893681 


£o,i[4,5] 


2.809797 


4.575161 


7.893681 




/ = 5 


1 = 10 


I = 20 


eo,i[2,2] 


10.981352 490 


17.965484 570 


29.949881 340088 


£o,j[2,3] 


10.981352 767 


17.965484 596 


29.949881 340530 


£ o,/[3, 3] 


10.981352 715 


17.965484 595 


29.949881 340562 


£o,j[3,4] 


10.981352 712 


17.965484 587 


29.949881 340559 


£o,i[4,4] 


10.981352 712 


17.965484 595 


29.949881 340559 


eo,i[4,5] 


10.981352 712 


17.965484 595 


29.949881 340559 
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Table 5: Comparison between PSLET, collecting the first 6 terms of the en- 
ergy series, £o,o[3, 3] Pade approximant and the results from direct numerical 
integration (DNI) [3]. 



a 


K = 6 


£ o,o[3, 3] 


DNI [3] 


1/100 


2.999906 259959 1 


2.999906 2599591 


- 


1/20 


2.997662 44644 


2.997662 44644 


2.997662 44644 


1/10 


2.990722 78231 


2.990722 78231 


2.990722 78232 


1/5 


2.963985 445 


2.963985 441 


2.963985 44193 


1/4 


2.944905 033 


2.944904 983 


2.944904 99229 


1/3 


2.906137 610 


2.906136 824 


2.906136 36892 


1/2 


2.809786 91 


2.809783 442 


2.809786 32134 


2 


1.932185 


1.933444 


1.932334 34201 
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